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Abstract: The quantum critical behavior of the Bose-Hubbard model for a description of two coupled 
Bose-Einstein condensates is studied within the framework of an algebraic theory. Energy levels, wave- 
function overlaps with those of the Rabi and Fock regimes, and the entanglement are calculated exactly 
as functions of the phase parameter and the number of bosons. The results show that the system goes 
though a phase transition and that the critical behavior is enhanced in the thermodynamic limit. 

PACS numbers: 73.43.Nq, 05.30.Jp, 03.75.Fi, 03.65.Ud 

As is well-known, the two-site Bose-Hubbard model^^' can be used to describe pair tunneling between 
two superconductors through an insulating junction, trapped ultra-cold bosonic gases, etc. Further- 
more, it can also be used to prepare macroscopically entangled states. The two-site Bose-Hubbard 
model has been investigated widely by many authors using various methods, such as the Gross-Pitaevskii 
approximation, ['^1 mean-field theory, the quantum phase model, and the Bethe ansatz method. 
In [8], the temporal evolution of the expectation value for the relative number of particles between the 
two condensates for different choices of the coupling parameter and distinct initial states is analyzed. 
Also, quantum phase transitions that occur at zero temperature as a function of a coupling constant have 
become important in connection with various quantum many-body systems, such as the quantum Ising 
and rotor models, Fermi liquids, '^"^l , and atomic nuclei. '^^1 There are distinct features in these systems 
at the critical points. The main purpose of the present paper is to study the critical behavior of the 
two-site Bose-Hubbard model as a function of the coupling parameter and the total number of bosons. 

Specifically, we consider the two-site Bose-Hubbard Hamiltonian with 

ii = -Ej{cU + Sc) + Ec{c^cc^c + dUdU), (1) 

where (c) and d^ (d) are boson creation (annihilation) operators in two traps or different hyperfine 
states. Ej is related to the Josephson coupling exchanging bosons between the two states, and Ec is 
related to the charging energy. In the present work we focus on the case where the effective interaction 
energy for the internal Josephson dynamics is negative, t^^l Ec < 0. In order to study the transitional 
patterns of the system, the Hamiltonian (1) is re-parameterized as 

4t 

H = H/Eq = -(1 - x)(cU + d'^c) [c^cc'^c + dUdU), (2) 

n -I- 1 

where Eq is a constant in arbitrary unit, n is the total number of bosons, and < a; < 1 is the phase 
parameter. The division of the second term in (2) by n -I- 1 serves to ensure that the boson rank of the 
two terms in H is the same, thereby making comparisons of results as a function of boson number more 
meaningful. For large values of n, this system is in the Rabi regime when a; 0, the Fock regime when 
a; ~ 1, and the Josephson regime when < x < 1. 

When X = 0, the system is in the Rabi regime with eigenstates given by 

\x = 0;n,,n,) = ^(cU dt)'H(,t _ dt)n.|o), (3) 

where n — rii + n2, and |0) is the boson vacuum state which is never degenerate. The corresponding 
eigenenergy of (2) is 

£'„i,„2(a; = 0) = n2 - ni. (4) 



When X = 1, the system is in the Fock regime with eigenstates given by 

^tni^tn2 

\x=l;ni,n2) = . — p^ |0), (5) 
which is two- fold degenerate under the permutation m #712- The corresponding eigenenergy is 

£„„„,(x = l) = -^(n2+n2). (6) 
In the Josephson regime, < a; < 1, the unitary transformation [13] 

c= -^{a-ib), d= -^{a + ib) (7) 
can be used to rewrite Hamiltonian (2) as 

H={1- x){b^b - ata) + —^S+{0)S-{0) - ^x, (8) 

where 

S+ (0) = b^^ + , 5- (0) = 52 + (9) 

are boson pairing operators. In a manner similar to what was done in [7] and [14], it can be shown that 
the Bethe ansatz eigenvectors used for diagonalizing the Hamiltonian (8) may be written as 

\x;C,n = 2k + i.,+ V2) = N S+ {y^P)S+ {y''P) ■ ■ ■ 1^2), (10) 

where Af is the normalization constant, |!^i,t'2) = a^'^^b^''^\0) with f j = or 1 for i = 1, 2 is the boson 
pairing vacuum state that satisfies 

a\i,U2) = b^\pi,U2) =0, (11) 

and 

S^iyP) = T^ + 7^' (12) 
1 - ^ 1 + yl''> 

in which {i = 1, 2, • • • , /c) arc spectral parameters that arc to be determined, and ( is an additional 
quantum number for distinguishing different eigenvectors with the same quantum number k. It can 
then be verified by using the corresponding eigen-equation that (10) is the solution when the spectral 
parameters j/-^^ (i = 1, 2, • • • , fc) satisfy the following set of equations: 

2x ( 2.2 + 1 2.1 + A 1-x ^ yf 

+ . , ..(C) ~ ..(0 ..(c) ..(c) ^^-^i 



for i= 1,2, - ■ ■ ,k with the corresponding eigen-energy given by 

Ei^\x) = (1 - ar) ^ 4) + (1 - - t^i) - ^x. (14) 

i—1 yi 

It should be noted that the solutions (13) and (14) are not valid in the a; = (Rabi) or 1 (Fock) limits. 

To explore transitional patterns, some low-lying energy levels of the system for x G [0, 1] with n = 10, 
40, 100, and 160, respectively, were calculated. The results, which arc shown in Fig. 1, clearly show that 
there is a minimum in the excitation energies around x ~ 0.25 — 0.35 which corresponds to the Josephson 
critical region. While for small boson numbers the critical region is not very well defined, it becomes 
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clearer and sharper with increasing n values. Also, the energy level density in the critical region increases 
with increasing n. 

To probe the nature of the critical point behavior more deeply, overlaps |(a;; n|a;o; n)| for the ground 
states with xq = or 1 for n = 10, 40, 80, and 120, respectively, were calculated. The results, which 
are given in Fig. 2, show that there is a crossover point at a certain nonzero amplitude for the overlaps 
|(a;; n\xo = 0; n)| and |(a;; n\xo = 1; when n is relatively small, which yields to a cross-over region with 
near zero amplitude when n becomes large. This critical region in the large n limit should be regarded as a 
two phase (Rabi and Fock) coexistence region in analogue to the situation occurring in other finite boson 
systems. '^^'^^1 Furthermore, there is a sharp change in |(x;n|a;o = 0;n)| around a critical point Xc ~ 0.25 
for large n. These results suggest that the largest absolute value of the derivative of |(a;;n|a;o = 0;n)| 
with respect to x occurs around the critical point in the thermodynamic limit. 

As is known, the entanglement measure for any pure bipartite system is defined by 

V = -Tr(Pc logjv Pc) = -Tr(pd log^v Pd), (15) 

where pc is the reduced density matrix obtained by taking the partial trace over the subsystem d, and 
similarly for pd, and N is the total number of the Fock states for given n. We use the logarithm to 
the base A'' instead of base 2 to ensure that the maximal measure is normalized to 1. Fig. 3 shows the 
entanglement measure of the system as a fimction of x for n = 10, 40, 80, and 120, respectively. It is 
obvious that there is always a peak in the measure at or near the critical point, which is consistent with 
the so-called critical point entanglement. I^^l For small n, the maximal value of the measure is near 1, while 
it decreases with increasing values of n. The peak grows sharper with increasing values of n, tracking the 
behavior of the overlap \ {x;n\xo = 0;n)\ and the excitation energies around the critical point. 

In summary, we have studied the quantum critical behavior of the two-site Bose-Hubbard model 
within the framework of an algebraic theory. Energy levels, overlaps of the wavcfunctions with those in 
both the Rabi and Fock regimes, and the entanglement were calculated exactly as a function of the phase 
parameter and the total number of bosons. The results show not only the quantum phase transition 
patterns of the model, but also that the critical behavior are greatly enhanced in the thermodynamic 
(large n) limit. This enhancement of critical phenomena with increasing of n should be common in 
other finite quantum many-body systems, as recently shown in the interacting boson model for atomic 
nuclei. [^''l More importantly, the results suggest that while the quantum phase transition of the two-site 
Bose-Hubbard model between the Rabi and the Fock regimes is rather smooth for small values of n, it 
becomes clearer and sharp in the thermodynamic limit. Therefore, such quantum phase transition should 
be observed macroscopically. 
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Fig. 2. Overlaps |{x;n|a;o;n)| for the ground states with xo = or 1 for n = 10, 40, 80, and 120, respectively. 
The full lines represent the curves of |(x;n|a;o = 0;n)|, and the dotted lines represent those of |(x;n|a;o = 
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Fig. 3. The entanglement measures of the model as a function of a; for n = 10, 40, 80, and 120, respectively. 
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